function [ blip, tout ] = get_blip_math_ts(amiga_model, ... 
    audio_filter_mode, filter_method, length)

filter_short_len = get_filter_short_len();
fs = get_sampling_fq();

[ts1, ts2] = get_amiga_filter_ts(amiga_model, audio_filter_mode);

if strcmpi(length, 'full') || strcmpi(filter_method, 'fq')
    tout = linspace(0, (fs-2)/fs, fs-1);
elseif strcmpi(length, 'short')
    tout = linspace(0, (filter_short_len-1) / fs, filter_short_len);
end

blip = impulse(ts1*ts2, tout)';
blip = fir_min_phase(blip);
irout = filter_impulse(amiga_model, blip, filter_method);

if strcmpi(length, 'short') && strcmpi(filter_method, 'fq')
    tout = tout(1:1:filter_short_len);
    blip = irout(1:1:filter_short_len);
end

blip = blip / trapz(blip);

end